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SUMMARY 


A  study  was  made  to  correlate  theoretical  stresses  and  strains  in  a  pressure_vessel 
having  a  reinforced  opening  with  experimental  data  obtained  for  progressive  stages  of 
elastic -plastic  deformation  up  to  the  point  of  failure.  The  pressure  vessel  consisted  of  a 
hemispherical  dome  12  inches  in  diameter  with  a  wall  thickness.;of  0.  060  inch  and  a  rein¬ 
forced  opening  1  inch  in  diameter  and  was  machined  from  6061  aluminum.  The  rein¬ 
forcement  was  an  integrally  machined  ring  1.  70  inches  in  outer  diameter  and  0.  300  inch 
thick.  In  general,  there  was  good  agreement^  between  the  experimental  data  and  the 
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theoretical  results  based  on  true  stress  -  cohventional  strain  data  and  modified  shell  di¬ 


mensions.  The  conclusion  was  drawn  that  the  theoretically  computed  elastic-plastic 
stresses  can  be  used  with  the  true  ultimate  tensile  strength  of  the  shell  material  to  give  a 


reasonably  accurate  prediction  of  the  burst  pressure.  Also,  the  elastic  stregs- 


concen^dtion  factor  can  be  used  together  with  the  true  ultimate  tensile  strength^),  give  ^ 
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a  conservative  approximation  of  this  burst  pressure 


INTRODUCTION 

Reinforcements  are  utilized  around  the  openings  in  shells  to  reduce  stresses  caused 
by  the  removal  of  material  from  the  parent  structure.  In  this  study,  experimental  and 
theoretical  analyses  of  the  elastic -plastic  stress  state  around  a  reinforced  opening  in  a 
spherical  pressure  vessel  were  obtainedfor  progressive  stages  of  deformation  up  to  failure 
The  most  common  measure  of  the  effectiveness  of  an  opening  reinforcement  is  an 
experimental  or  theoretical  elastic  stress-concentration  factor.  Usually  this  factor  is 
combined  with  the  yield  strength  of  the  shell  material  to  determine  the  pressure  load  at 
which  incipient  local  plastic  flow  takes  place.  |~An  ideal  reinforcement  is  one  in  which  the 


pressure  at  which  local  yielding  (plastic  flow)  occurs  is  the  same  as  if  the  pressure  ves¬ 
sel  did  not  contain  an  opening  or  reinforcement.] 

Elastic  stress-concentration  factors  have  been  investigated  for  apex  openings  in  bi- 
axially  stressed  plates  (refs.  1  to  5).  /  A  small  amount  of  local  plastic  flow  can  very  sig¬ 
nificantly  reduce  the  stress-concentration  factor] (ref.  6^;  therefore,  the  elastic  stress- 
concentration  factor  may  result  in  an  unduly  conservative  prediction  of  the  ability  of  ves¬ 
sels  with  reinforced  openings  to  withstand  internal-pressure  loadings.  It  therefore  ap¬ 
pears  that  an  elastic -plastic  stress  analysis  is  required  to  provide  a  more  realistic  eval¬ 
uation  of  the  effectiveness  of  the  reinforcement J  , 

References  7  and  8  present  numerical  analyses  of  a  continuous  elastic-plastic  shell 
that  can  be  modified  to  treat  the  problem  of  a  reinforced  opening.  Reference  9  presents 
a  numerical  analysis  that  is  directly  applicable  to  the  reinforced  opening  problem.  All 
of  these  theoretical  methods  are  based  on  the  von  Mises  yield  criterion  and  associated 
flow  rules  and  on  the  assumption  of  infinitesimal  strains.  In  addition,  shear  deforma¬ 
tions  are  neglected.  The  validity  of  these  assumptions  requires  experimental  verifica¬ 
tion  for  the  case  of  the  finite  deformation  of  a  shell  that  contains  a  discontinuity  such  as 
an  opening  reinforcement. 

The  purposes  of  this  investigation  are  threefold:  first,  to  evaluate  the  accuracy  with 
which  plasticity  theory  and  a  numerical  shell  analysis  can  predict  the  elastic -plastic 
stress  and  strain  states  in  a  discontinuous  shell  structure  loaded  to  failure;  second,  to 
evaluate  the  validity  of  the  elastic  stress -concentration  factor  as  a  measure  of  the  ef¬ 
fectiveness  of  a  reinforcement  around  an  opening  in  a  pressure  vessel;  third,  to  use 
these  results  to  draw  some  preliminary  conclusions  about  a  method  for  predicting  the 
burst  pressure  and  the  mode  of  failure  for  such  a  reinforced  vessel. 

Experimental  strain  data  to  the  burst  point  were  obtained  from  a  hemispherical  dome 
with  a  nominal  diameter  of  12  inches,  a  wall  thickness  of  0.  060  inch,  and  a  reinforced 
circular  apex  opening  1  inch  in  diameter.  The  reinforcement  proportions  were  those  of 
the  most  effective  configuration  for  strengthening  a  hole  in  a  biaxially  stressed  sheet  with 
a  stress  field  ratio  of  1,  as  reported  in  reference  2.  The  strain  data  were  correlated 
with  results  from  a  conventional  elastic-edge-influence-coefficient  analysis  and  the 
elastic-plastic  shell  analysis  of  reference  9  and  are  presented  in  this  report. 


SYMBOLS 


a 

radius  of  dome,  in. 

h,  hT 

thickness  of  shell  and  ring, 

E 

modulus  of  elasticity,  psi 

respectively 

H,V 

radial  and  axial  stress  re¬ 

M 

couple,  lb-in.  /in. 

sultants,  lb/in. 

P 

pressure,  lb/sq  in. 
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r 

radial  coordinate,  in. 

r 

radial  coordinate  to  centroid 
of  ring,  in. 

u 

radial  displacement,  in. 

P 

meridional  rotation,  rad 

influence  coefficients  for  de¬ 
flection  and  rotation,  re¬ 
spectively 

e 

strain,  in.  /in. 

V 

Poisson’s  ratio 

<7 

stress,  psi 

a 

true  stress,  psi 

ae 

effective  stress, 

V«  +  <Te’  pS‘ 

meridional  slope,  rad 

Subscripts: 

e 

effective 

et,  ep 

effective  total,  effective 
plastic 

H,  M 

caused  by  a  unit  radial 
force  and  meridional 
couple,  respectively 

0 

original 

max 

maximum 

mem 

membrane 

P 

caused  by  a  unit  internal 
pressure 

R,  S 

ring,  shell 

t 

transverse  direction 

u 

uniaxial  test  data 

y 

yield 

£,£P 

transverse,  transverse 
plastic 

0,  Be,  flp 

circumferential,  circum¬ 
ferential  elastic,  cir¬ 
cumferential  plastic 

cp,(pp 

meridional,  meridional 
plastic 

1,2 

inner  and  outer  circum¬ 
ference  of  reinforcement, 
respectively 

EXPERIMENTAL  APPARATUS  AND  PROCEDURE 
Test  Specimen 

The  geometry  of  the  hemispherical-dome  test  specimen  is  shown  in  figure  1.  The 
dome  had  a  nominal  diameter  of  12  inches,  a  wall  thickness  of  0. 060  inch,  and  an  apex 
opening  1  inch  in  diameter.  The  reinforcement  ring  was  integrally  machined  to  an  outer 
diameter  of  1.  70  inches,  a  thickness  of  0.  30  inch,  and  had  a  0.  030-inch  radius  fillet  at 
the  junction  with  the  shell.  Fabrication  tolerances  were  specified  as  ±0.  002  inch  for  the 
dome  and  ring  diameters  and  ±0.  001  inch  for  the  dome  wall  and  ring  thicknesses.  The 
dome  and  four  tensile  coupons  were  cut  from  a  billet  of  6061-0  aluminum  alloy  and  heat 
treated  to  the  T6  condition  prior  to  final  machining. 
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Figure  2(a)  shows  the  method  of  closing  the  base  of  the  dome.  A  removable  head  was 
attached  and  sealed  with  an  alloy  that  melted  at  about  300°  F,  as  described  in  refer¬ 
ence  10.  A  cooling  coil  was  wrapped  around  the  dome  to  protect  the  strain  gages  from 
the  heat  during  installation  of  the  removable  head.  The  apex  opening  of  the  dome  was 
sealed  by  a  stainless  steel  plug  and  a  rubber  O-ring,  as  shown  in  figure  2(b). 

Strain-Gage  Installation 

Two  types  of  electrical  resistance  strain  gages  were  used  in  this  investigation. 

First,  conventional  Advance  foil  gages  with  l/8-inch-square  active  elements  were  in- 
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stalled  to  measure  strains  below  the  elastic  limit.  These  were  replaced  later  by  post¬ 
yield  Advance  gages  of  the  same  size  to  measure  strains  in  the  elastic -plastic  range. 

The  postyield  gages  that  were  used  were  found  from  tests  to  be  less  accurate  than  the 
conventional  gages  in  the  elastic  range  and  were,  therefore,  not  used  to  measure  the 
elastic  strains  up  to  the  yield  point. 

Twenty-eight  conventional  gages  and  28  postyield  gages  were  mounted  at  the  locations 
indicated  in  figure  1(b).  Two  circumferentially  oriented  gages  and  two  axially  oriented 
gages  were  placed  along  the  inner  circumference  of  the  reinforcing  ring  at  location  A. 
Three  pairs  of  gages  were  installed  on  both  the  inside  and  outside  surfaces  of  the  shell 
with  the  gage  centerlines  0.  1  inch  from  the  reinforcement  ring  (location  B).  Each  pair 
consisted  of  meridionally  and  circumferentially  oriented  gages.  Installations  similar  to 
location  B  were  made  at  location  C,  approximately  midway  between  the  reinforcement 
and  the  end  closure.  Gages  were  bonded  to  the  specimen  with  cyanoacrylate  cement. 

Lead  wires  from  the  interior  gages  were  brought  through  the  tank  base  by  pass-through 
connectors  (fig.  2(a)).  Lead  wires  from  location  A  passed  through  the  split  ring  in  the 
closure  plug  assembly  as  shown  in  figure  2(b).  A  multichannel  digital  recorder  utilizing 
a  five-wire  system  was  used  to  measure  and  tabulate  the  gage  output. 


Testing  Procedure 

The  dome  was  loaded  by  internal  hydraulic  pressure  that  was  measured  on  Bourdon 
gages  having  a  maximum  error  of  0.  5  pound  per  square  inch.  In  the  elastic  range,  gage 
readings  were  taken  at  pressures  of  0,  100,  200,  300,  400,  500,  and  600  pounds  per 
square  inch.  Each  pressure  increment  resulted  in  a  membrane  stress  increase  of  about 
5000  pounds  per  square  inch.  The  pressure  was  released  after  each  reading,  and  a  zero 
reading  was  again  taken.  A  uniform  time  interval  of  4  minutes  was  maintained  between 
readings  to  minimize  viscoelastic  effects,  and  a  duplicate  run  was  made.  The  maximum 
difference  between  the  two  runs  was  less  than  1  percent. 

Following  completion  of  testing  in  the  elastic  range,  the  dome  was  reinstrumented 
with  postyield  gages  and  pressurized  to  burst.  Data  were  taken  in  the  course  of  this 
loading  at  pressures  of  0,  600,  650,  700,  750,  800,  850,  and  900  pounds  per  square  inch. 
The  load  was  held  constant  at  each  of  these  values  only  long  enough  to  record  readings 
from  all  28  gages,  in  an  effort  to  minimize  the  effects  of  creep,  which  were  far  more  ser¬ 
ious  than  viscoelastic  effects  in  the  plastic  range.  This  time  interval  was  somewhat 
under  a  minute.  Bursting  of  the  specimen  occurred  at  a  pressure  of  935  pounds  per 
square  inch  about  10  minutes  after  the  postyield  loading  was  begun. 
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Reduction  of  Data 


Data  from  the  conventional  gages  were  entered  into  a  computer  program  for  the  re¬ 
duction  of  elastic -strain  readings,  which  corrected  them  for  zero  drift  and  nonlinearity  by 
the  method  of  least  mean  squares  and  converted  the  readings  to  stresses  and  strains. 

The  data  from  the  two  runs  were  averaged.  Output  data  from  the  postyield  gages  were 
processed  by  a  second  program  that  used  an  iterative  procedure  to  calculate  the  elastic 
and  plastic  components  of  the  measured  strains.  By  using  a  stress-strain  curve  based 
on  the  true  stress  and  conventional  strain  test  data  obtained  from  the  tensile  coupons 
machined  from  the  billet,  stresses  in  the  plastic  range  were  calculated  as  explained  in 
appendix  A. 


THEORETICAL  ANALYSIS 
Elastic  State 


Stresses  in  the  shell  prior  to  yielding  were  calculated  by  a  standard  edge-influence- 
coefficient  technique  and  also  by  a  numerical  solution  of  the  governing  differential  equa¬ 
tions.  The  first  method  is  explained  in  detail  in  appendix  B,  and  the  second  was  taken 
from  reference  9. 

Influence-coefficient  analysis.  -  A  number  of  methods  are  available  for  calculating 
the  edge  deformations  in  a  shallow  shell  as  a  result  of  pressure  and  boundary  loads. 
Galletly  (ref.  4)  analyzed  three  of  these  methods  and  recommended  Esslinger's  solution, 
in  which  the  approximation  cot  cp  =  l ftp  is  made.  This  approximation  is  acceptable  for 
values  of  the  meridional  slope  cp  up  to  about  30°,  for  which  the  error  would  be  about 
10  percent.  Because  of  the  local  nature  of  the  effects  of  boundary  forces  and  moments 
in  shells,  significant  values  of  cp  would  be  limited  to  about  30°  in  a  dome  with  an  apex 
opening  half-angle  of  about  15°.  In  the  present  investigation,  the  spherical  dome  con¬ 
tained  an  apex  opening  with  a  half-angle  of  about  8°  as  measured  to  the  outside  diameter 
of  the  reinforcement.  Therefore,  Esslinger's  coefficients  were  applicable.  Combined 
with  the  well-known  ring-deformation  equations  (ref.  11),  they  yielded  discontinuity  for¬ 
ces  and  moments  at  the  ring-shell  junction,  which  in  turn  led  to  the  elastic  stresses. 

Numerical  analysis.  -  In  an  alternative  elastic-state  analysis  based  on  reference  9, 
the  reinforcement  was  treated  as  a  short  spherical  segment.  Equations  of  equilibrium 
and  compatibility  for  both  the  reinforcement  and  the  dome  segments  were  reduced  to  two 
second-order  ordinary  differential  equations.  The  independent  variable  in  these  equa¬ 
tions  was  the  meridional  slope  cp.  The  two  dependent  variables  were  the  change  in  mer¬ 
idional  slope  and  the  radial  stress  resultant.  These  differential  equations  were  based  on 
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the  usual  thin-shell  assumptions,  namely,  that  transverse  stresses  and  shearing  strains 
are  negligible  and  that  wall  thicknesses  are  much  smaller  than  radii  of  curvature. 

Boundary  conditions  on  the  basic  differential  equations  were  as  follows.  At  the 
inner  circumference  of  the  reinforcement,  the  radial  stress  resultant  vanished,  while 
the  axial  stress  resultant  and  the  meridional  couple  together  were  statically  equivalent 
to  the  closure-plug  load  as  applied  at  the  centerline  of  the  0-ring  a  small  distance  from 
the  edge  (fig.  2(b)).  At  the  common  boundary  between  the  reinforcement  segment  and  the 
shell  segment,  corresponding  edge  forces  and  deformations  in  the  two  segments  were 
made  equal.  Finally,  at  the  outer  boundary  of  the  dome,  which  was  assumed  to  lie  in  a 
membrane  region,  the  transverse  stress  resultant  and  the  meridional  couple  vanished. 

The  reinforcement  and  the  dome  meridians  were  subdivided  into  equal  increments 
0.  05  inch  in  length.  The  two  differential  equations  were  expressed  in  finite -difference 
form,  and  a  numerical  integration  technique,  which  is  given  in  detail  in  reference  9,  was 
applied  to  solve  them.  The  meridional-slope  changes  and  the  radial  stress  resultants 
thus  obtained  were  used  to  determine  the  elastic  stress  and  strain  states  throughout  the 
shell  and  reinforcing  ring. 


Elastic-Plastic  State 

The  numerical  analysis  of  reference  9  was  also  applied  to  the  shell  in  the  elastic- 
plastic  state.  As  explained  in  appendix  A,  the  deformation  theory  of  plasticity  and  the 
method  of  successive  approximations  were  used  to  determine  the  plastic  components  of 
strain.  The  solution  of  the  basic  differential  equations  for  each  approximation  was  iden¬ 
tical  to  the  elastic -state  solution.  Boundary  and  junction  conditions  were  modified  to  in¬ 
clude  the  effects  of  plastic  strain  and  were  then  incorporated  into  the  numerical  solution, 
as  explained  previously. 

The  elastic-plastic  analysis  was  applied  to  the  reinforced  dome  under  conditions 
ranging  from  incipient  plastic  flow  to  burst.  It  was  therefore  necessary  to  make  some 
provision  for  the  small  but  finite  middle  surface  strains,  which  would  certainly  be  pres¬ 
ent  under  this  wide  range  of  loads.  This  was  done  by  assuming  that,  at  a  given  pressure, 
the  curvature  and  the  wall  thickness  throughout  the  shell  were  equal  to  the  curvature  and 
the  wall  thickness  in  the  membrane  region  of  the  dome.  These  latter  dimensions  are 
calculated  in  appendix  C  (see  table  I,  p.  21).  The  true  effective  stress  -  conventional 
effective  strain  curve  that  was  based  on  the  data  in  table  I  is  shown  in  figure  3  together 
with  an  exponential  curve  fit.  This  use  of  the  instantaneous  membrane  dimensions  to¬ 
gether  with  true  stress  -  conventional  strain  data  ensured  that  the  correct  theoretical 
stress  -  strain  -  pressure  relations  would  exist  in  the  membrane  region  of  the  dome  even 
with  small  finite  strains. 
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Figure  3.  -  Effective  stress  -  effective  total  strain  relations. 


RESULTS  AND  DISCUSSION 
Elastic  State 

The  results  of  the  experimental  and  theoretical  elastic-stress  analyses  of  the  dome 
test  specimen  are  compared  in  figure  4.  Stresses  are  normalized  with  respect  to  the 
nominal  membrane  stress  in  the  shell.  Circumferential  and  meridional  stresses  that 
were  calculated  by  the  numerical  method  of  reference  9  are  presented  for  both  inside  and 
outside  surfaces  along  a  meridian.  Experimental  data  are  represented  by  rectangular 
symbols  whose  dimensions  indicate  the  gage  length  and  the  scatter.  Agreement  between 
the  experimental  data  and  the  theoretical  stresses  is  generally  good. 

The  maximum  theoretical  stress  occurred  on  the  inside  surface  of  the  shell  at  the 
ring-shell  junction.  At  this  location  the  nondimensional  meridional  stress  was  1.  42  from 
calculations  based  on  either  reference  9  or  the  influence-coefficient  analysis,  which  is 
explained  in  appendix  B.  This  stress  concentration  was  caused  mainly  by  the  discontin¬ 
uity  moment  at  the  junction  because  the  effective  stress-concentration  factor  of  the  mid¬ 
dle  surface  was  0.  975. 
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Figure  4.  -  Experimental  and  theoretical  elastic  stress. 

Elastic-Plastic  State 

Curves  of  strain  in  the  dome  as  a  function  of  internal  pressure  are  shown  in  figure  5 
for  locations  B  (near  the  ring-shell  junction)  and  C  (in  the  membrane  region).  The  ex¬ 
perimental  data  are  shown  for  pressures  above  600  pounds  per  square  inch,  the  pressure 
for  incipient  plastic  flow,  and  terminate  at  the  last  pressure  at  which  data  were  taken  be¬ 
fore  the  strains  went  out  of  the  range  of  the  measuring  equipment.  The  experimental 
burst  pressure  is  also  indicated.  Most  of  the  experimental  and  theoretical  data  become 
increasingly  nonlinear  at  pressures  above  800  pounds  per  square  inch  and  increase  so 
rapidly  as  to  indicate  an  instability  condition  above  900  pounds  per  square  inch.  For  both 
the  inside  surface  (fig.  5(a))  and  the  outside  surface  (fig.  5(b))  at  location  B  near  the 
reinforcement,  the  theoretical  strains  are  substantially  lower  than  the  measured  strains. 
The  greatest  discrepancy  between  experiment  and  theory  is  at  the  outside  surface  at  lo¬ 
cation  B,  as  is  shown  in  figure  5(b)  by  the  meridional  strain  curves.  The  theoretical 
analysis  predicted  a  strain  that  decreased  at  pressures  above  800  pounds  per  square  inch 
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Strain,  in. /in. 


.024 


(a)  Location  B,  inside  surface. 


Pressure,  Ib/sqin. 


(c)  Location  C  (membrane  region),  inside  and  outside  surfaces. 

Figure  5.  -  Experimental  and  theoretical  strains  in  elastic-plastic  state.. 


and  remained  in  the  elastic  state, 
while  the  measured  strains  actually 
increased  and  experienced  substantial 
plastic  flow.  In  the  membrane  region 
(fig.  5(c))  the  wide  scatter  band  in¬ 
cludes  data  from  12  gages  oriented  in 
meridional  and  circumferential  direc¬ 
tions  on  both  inside  and  outside  sur¬ 
faces  in  the  nominal  membrane  region 
of  the  shell.  This  scatter  was  prob¬ 
ably  caused  mainly  by  distortion  of 
the  dome  and  nonuniform  reduction  in 
wall  thickness  with  increasing  pres¬ 
sure. 

In  figure  6  experimental  and  theo¬ 
retical  results  are  correlated  on  the 
basis  of  true  stress  as  a  function  of 
pressure.  Because  stresses  increase 
slowly  with  increasing  strain  in  the 
plastic  state,  agreement  between  ex¬ 
perimental  and  theoretical  results  is 
better  for  stresses  than  for  strains. 
The  experimental  burst  pressure  and 
the  theoretical  membrane  burst  pres¬ 
sure  from  table  I  (p.  21)  are  also  in¬ 
dicated. 

Both  the  experimental  and  the 
theoretical  data  presented  in  figure  6 
show  that  the  material  at  location  A 
remains  elastic  up  to  the  burst  pres¬ 
sure.  The  theoretical  solution  shows 
a  tensile  instability  in  the  membrane 
region  at  a  pressure  of  969  pounds 
per  square  inch.  The  experimental 
data  for  location  C  show  that  the 
stress  field  ratio  of  1  was  not  main¬ 
tained  as  plastic  flow  increased.  At  a 
pressure  of  900  pounds  per  square 
inch  the  circumferential  stresses  were 
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Figure  6.  -  Experimental  and  theoretical  stresses  in  elastic-plastic  state. 


generally  higher  than  the  meridional  stresses.  The  maximum  measured  stresses  on  the 
inside  surface  at  location  B  in  the  meridional  direction  were  higher  than  those  predicted 
for  this  location.  The  theoretical  curve  for  this  location  reached  the  ultimate  tensile 
strength  of  the  material  at  a  pressure  of  933  pounds  per  square  inch,  however,  or  only 
0.  2  percent  below  the  experimental  burst  pressure  of  935  pounds  per  square  inch. 


A1  0  Burst  Condition 

j;.|  t  YSJ^.  ‘ 

A  view  of  the  dome:  fracture  and  meridional  sections  through  the  junction  region  are 
shown  in  figure  7.  pThe^failure  apparently  originated  on  the  inside  surface  approximately 
0.  1  inch  from  the  ring- shell  junction.)  The  rough  surface  of  the  dome  apparent  in  fig- 
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ure  7(a)  was  caused  by  high  plastic  flow. 
The  specimen  had  an  initial  root  mean 
square  surface  finish  of  32  microinches. 

The  sectional  views  in  figure  7(b) 
show  a  definite  necking  in  the  dome  ma¬ 
terial  a  small  distance  from  the  ring- 
shell  junction.  There  is  no  apparent 
thinning  at  the  junction  itself,  however, 
because  of  the  transverse  restraint  of  the 
ring  and  fillet.  This  restraint,  together 
with  the  slight  thickening  of  the  dome  wall 
near  the  root  of  the  fillet,  probably 
caused  the  maximum  stress  to  be  dis¬ 
placed  from  the  fillet  to  the  point  at  which 
the  shell  wall  assumed  a  uniform  thick¬ 
ness. 

The  experimental  burst  results  and 
the  theoretical  maximum  stresses  are 
correlated  in  figure  8.  The  stresses  are 
nondimensionalized  with  respect  to  the 
nominal  membrane  stress  for  the  dome 
computed  from  the  original  shell  dimen¬ 
sions.  If  no  fillet  had  been  present,  the 
highest  stress  would  have  occurred  at 
the  ring-shell  junction  on  the  inside  sur¬ 
face  of  the  shell.  When  a  small  fillet  is 
present,  it  is  reasoifable  to  assume  that 
the  critical  location  is  transferred  to  the 
end  of  the  fillet. 

)  The  elastic  stress- concentration 

j 

factor  for  the  end  of  the  fillet  on  the  in¬ 
side  surface  is  1.  36  from  the  analysis  of 
reference  9.  On  the  basis  of  the  ultimate 
tensile  strength,  this  would  lead  to  a  pre¬ 
diction  of  a  burst  pressure  of  813  pounds 
per  square  inch,  which  is  13  percent  be¬ 
low  the  experimental  burst  pressure,]  ^ 
The  nondimensional  el  as  tic -pi  as  tic  stress 
curve  from  the  analysis  of  reference  9 
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Figure  8.  -  Correlation  of  experimental  burst  results  and  theoretical  predictions  for  burst  pressure. 
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gives  a  predicted  burst  pressure  of  894  pounds  per  square  inch,  about  4  percent  below 
the  experimental  burst  pressure.  The  predicted  burst  pressure  for  the  membrane  region 
of  the  dome  in  which  failure  would  occur  through  tensile  instability  (table  I,  p.  21)  is 
969  pounds  per  square  inch,  about  4  percent  above  the  experimental  burst  pressure. 

{ Based  on  the  maximum  stress  failure  theory,  the  elastic  stress -concentration  fac¬ 
tor  gives  the  minimum  predicted  burst  strength  of  the  shell,  while  the  elastic-plastic 
membrane  instability  analysis  gives  the  theoretical  maximum  pressure  the  shell  can  with¬ 
stand^  The  difference  between  these  two  values  would  determine  whether  it  would  be 
worthwhile  to  go  through  an  elastic -plastic  shell  analysis.  The  membrane  stress  curve 
for  a  shell  without  a  hole  intersects  the  ultimate  tensile  stress  curve  at  a  nondimensional 
stress  of  1.  14.  This  is  the  elastic  stress-concentration  factor  that  the  dome  could  toler¬ 
ate  without  any  deleterious  effect. 


i  conclusions/ 

The  following  general  conclusions  are  drawn  from  the  results  of  this  investigation  of 
the  elastic -plastic  stress  and  strain  states  in  a  pressurized  spherical  shell  with  an  apex 
opening  reinforced  by  an  integrally  machined  ring: 

'PT?  The  numerical  shell  analysis  and  the  shallow-shell  edge-influence-coefficient 
analysis  give  identical  results  for  the  elastic  stress -concentration  factor  in  a  shell  con¬ 
taining  a  discontinuity  such  as  a  reinforced  opening. 

2.  Correlation  between  experimental  and  theoretical  stresses  and  strains  in  the 
elastic  state  is  very  good. 
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I  3.  The  deformation  theory  of  plasticity  and  the  method  of  successive  approximations 
can  be  combined  with  a  numerical  shell  analysis  to  calculate  elastic-plastic  stresses 
that  are  generally  in  good  agreement  with  experimental  results. 

4.  The  elastic-plastic  stresses  based  on  true  stress  data  and  the  instantaneous  mem¬ 
brane  shell  dimensions  could  be  used  with  the  true  ultimate  tensile  strength  of  the  shell 
material  to  give  a  reasonably  accurate  prediction  of  the  burst  pressure  in  a  shell  con¬ 
taining  a  discontinuity. 

5.  The  elastic  stress-concentration  factor  may  be  used  in  conjunction  with  the  true 
ultimate  tensile  stress  for  a  ductile  material  to  give  a  conservative  approximation  of  the 
burst  pressure.  The  discrepancy  between  the  experimental  burst  pressure  and  the 
burst  pressure  calculated  from  the  elastic  stress -concentration  factor  was  13  percent 

in  this  investigation. 

.....  ^  ,/ 

Lewis  Research  Center, 

National  Aeronautics  and  Space  Administration, 

Cleveland,  Ohio,  December  10,  1964. 


APPENDIX  A 


ELASTIC-PLASTIC  STRESS  AND  STRAIN  EQUATIONS 

Derivations  of  equations  similar  to  or  identical  with  the  following  equations  are  pre¬ 
sented  in  detail  in  reference  12  and  will  not  be  repeated  here.  The  stress  and  strain 
equations  are 


-.yH 


V 


eet =  ~3~  \(<>-ee>2  +  (ee  ~h)2  +  {erecp)2 


(Ala) 


_  % 

e  ep  =  e  u  g 


et 


for  a0  <  ay 

_2il±u)^e  for  > 

3  e  e  y 


(Alb) 


in  which  the  subscript  u  indicates  the  results  of  a  uniaxial  tensile  test,  and  the  effective 
stress  ae  is  a  known  function  of  as  determined  by  a  curve  fit  of  experimental  data 
(fig.  3).  The  plastic  components  of  strain  are  then  given  as 


e<pp 


-  e 


(A2a) 


e0P  3eeP  (2e0  ‘  U  ‘  ecp) 
et 


Up~  ~(e<pp  +  e&P} 


(A2b) 

(A2c) 


The  total  transverse  strain  is 

er-Th^[u(ecP+e0)  +  (1'2l')(e<pp  +  e0P)]  (A3) 


Finally 


cr 


<P 


■v 


+  v(e 


(A4a) 
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The  successive  approximation  technique  for  reducing  experimental  strain  data  in  the 
elastic-plastic  state  is  as  follows.  The  measured  total  strains  e  and  eQ  are  com¬ 
bined  with  an  assumed  total  transverse  strain  e  ^  in  equation  (Ala)  to  produce  a  trial 
value  of  the  equivalent  total  strain  e^.  By  using  this  trial  value,  first  approximations 
to  the  equivalent  plastic  strain  e  and  the  individual  plastic  strains  e  ^  and  €  can 
be  calculated  by  using  equations  (Alb)  and  (A2).  A  second  approximation  to  the  trans¬ 
verse  strain  e  ^  is  then  calculated  from  equation  (A3),  and  the  procedure  is  repeated 
until  the  desired  degree  of  convergence  is  obtained.  Stresses  are  then  calculated  by 
using  equations  (A4). 

A  similar  procedure  is  used  in  the  shell  analysis  of  reference  9.  In  this  analysis, 
the  basic  differential  equations  that  describe  the  conditions  of  equilibrium  and  compata- 
bility  of  deformation  throughout  the  shell  contain  the  plastic  strain  terms  e  ^  and  e  ^p. 
For  a  first  approximation,  these  terms  are  set  equal  to  zero  throughout  the  shell*  Solu¬ 
tion  of  the  basic  differential  equations  then  yields  the  elastic -state  strains.  These 
strains  are  substituted  into  equation  (Ala),  and  the  resulting  value  of  is  used  to  cal¬ 
culate  €  ep  and  a  second  approximation  to  the  plastic  strains  e  ^  and  e  ^ .  The  basic 
differential  equations  are  re-solved  with  these  improved  plastic  flow  terms,  and  the 
process  is  repeated  until  convergence  is  reached.  Stresses  are  then  computed  by  means 
of  equations  (A4). 
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APPENDIX  B 


INFLUENCE-COEFFICIENT  ANALYSIS 

As  shown  by  figure  9,  compatibility  of  deformations  at  the  junction  of  the  ring  and 
shell  requires  that 


In  influence-coefficient  notation,  equations  (Bla)  become 


(Bla) 


R,  H  l  2  C°S  <f>2  +  H2 )  +  §R,  MMR  “  6S,  HH2  +  6S,  MM2  +  6S,  pp 


V  H  (  2  cos  ^2  +  H2 )  +  eR,  MMR  "  6S,  HH2  +  0S,  MM2  +  dS,  pP 


(Bib) 


Figure  9.  -  Force  diagram  for  influence-coefficient  analysis. 


The  closure  plug  load  and  the  axial 
stress  resultant  Vg  are 


(B2) 

The  resultant  twisting  moment  on  the  ring,  re¬ 
ferred  to  the  centroid  of  its  cross  section,  is  then 

MR  =  Y  [r2M2  '  (rl  +  Arl)(F  -  rl  -  Arl>Vl 

-r2^r2  "  F)V2] 

For  the  dimensions  shown  in  figure  1 
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Mr  =  1.2593  M2  -  0.  1145  p 


(B3) 


From  reference  11,  the  ring  influence  coefficients  for  E  =  10*  pounds  per  square 
inch  and  v  =  0.  333  inch  per  inch  are 


2  r2 

5  =  —  (  _? - -  -  v\  =  4.  887x10' 7  in.  /(lb/in. ) 

R’ H  Eh’  \  2  2 

r2  '  rl 


.  = - — -  =  565.  4xl0~7  rad/((lb-in.  )/in. ) 

JX)  jVI 


E(h')3  InL2 


>  (B4a) 


0R,H  "  6R,M  0R,M 


fa2  -  r2  -  |42  -  r2)  =  15.  66xl0"7  rad/(lb/in. ) 


or  in.  /((lb-in.  )/in. ) 

From  reference  5,  the  shell  influence  coefficients  for  <p2  =  0.  1414  radian  are 

6g  H  =  -59.  909X10"7  in.  /(lb/in. ) 

6g  M  =  652.  79X10' 7  in.  /((lb-in.  )/in. ) 

5C  „  =  28.489X10"7  in. /psi 

D,  P 

0g  H  =  652.  79X10"7  rad/(lb/in. ) 

0g  m  =  -20  321X10'7  rad/((lb-in.  )/in. ) 


Vs0 

Substitution  of  equations  (B3)  and  (B4)  into  (Bib)  yields  the  following  discontinuity 
force  and  couple: 


(B4b) 


Hx  =  0.  3339  p 


•\ 


Mj  =  0.  009353 


. 


(B5) 
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Finally,  the  elastic  stress- concentration  factor  for  the  dome  is 


pa  |  H2  cos  ^2  | 


=  1. 420  p 
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APPENDIX  C 


ELASTIC-PLASTIC  MEMBRANE  ANALYSIS 

The  first  step  in  the  calculation  of  stresses,  strains,  and  dimensions  in  the  mem¬ 
brane  region  of  the  spherical  dome  (table  I)  is  the  conversion  of  uniaxial  tensile  test  data 
to  true  effective  stress  as  a  function  of  conventional  strain.  If  A  and  Aq  represent  the 
instantaneous  and  original  cross-sectional  areas,  respectively,  of  the  tensile  specimen, 
then 

A  =  A0(l+et)2 

in  which  the  subscript  t  denotes  the  transverse  direction.  If  the  material  is  assumed  to 
be  incompressible  in  both  the  elastic  and  the  plastic  states,  the  true  uniaxial  stress  be¬ 
comes 


a 

u 


(1  -  *u/2>2 


TABLE  I.  -  CALCULATION  OF  MEMBRANE  DIMENSIONS  AND  BURST  PRESSURE 


Uniaxial  tensile  data 

True 

Effective 

Effective 

Elastic 

Plastic 

Total 

Total 

Dome 

radius 

Dome 

wall 

Internal 

pressure, 

Conven¬ 

tional 

stress, 

V 

psi 

Conven¬ 

tional 

strain, 

eu’ 

in.  /in. 

effective 

stress, 

V 

psi 

plastic 

strain, 

eep’ 
in.  /in. 

total 

strain, 

eet’ 
in.  /in. 

circum¬ 

ferential 

strain, 

6ee’ 
in.  /in. 

ferential 

strain, 

eBp’ 
in.  /in. 

ferential 

strain, 

ee> 

in.  /in. 

LX  CUiD  V  vl  uL 

strain, 

in.  /in. 

a, 

in. 

thickness, 

h, 

in. 

P, 

lb/sq  in. 

36  000 

0.00360 

36  100 

0 

0.  00321 

0.00241 

0 

0. 00241 

-0.  00241 

6.  044 

0.  0599 

716 

38  000 

. 00395 

38  200 

.00013 

. 00353 

.00255 

. 00006 

.00261 

-.00268 

6.046 

.0598 

756 

40  000 

. 00470 

40  200 

.00068 

. 00425 

. 00268 

. 00034 

. 00302 

-.00336 

6.048 

.0598 

795 

42  000 

. 00595 

42  300 

.00173 

. 00549 

.00281 

. 00086 

. 00367 

00454 

6.052 

.0597 

835 

44  000 

. 00830 

44  400 

.00386 

.00781 

.00296 

.00193 

. 00489 

-.00682 

6.059 

.0596 

873 

46  000 

.0172 

46  800 

.0125 

.0167 

.00312 

. 00625 

. 00937 

-.0156 

6. 087 

.0591 

909 

48  000 

.0323 

49  600 

.0273 

.0317 

.00331 

.  0136 

.0169 

-.0306 

6.132 

.0582 

941 

49  000 

.0429 

51  100 

.0378 

.0423 

.00341 

.0189 

.0223 

-.0412 

6. 164 

.0575 

953 

50  000 

.0560 

52  800 

.0507 

.0554 

. 00352 

.0254 

.0289 

-.  0542 

6.  204 

.0567 

965 

50  300 

.0614 

53  400 

.0561 

.0609 

.00356 

.0280 

.0316 

-.0597 

6.  220 

.0564 

968 

50  500 

.0668 

53  900 

.0614 

.0662 

. 00359 

.0307 

.0343 

-.0650 

6.237 

.0561 

a969 

50  700 

.0770 

54  600 

.0715 

.0764 

. 00364 

.0358 

.  0394 

-.0759 

6.  268 

.0555 

967 

aMembrane  burst  pressure. 


Since  equation  (A4c)  shows  that  the  effective  stress  is  equal  to  the  nonzero  stress  in  a 
uniaxial  stress  state,  with  second-order  strain  terms  neglected,  the  true  effective  stress 
becomes 


a  =  a  (1  +  e  )  (Cl) 

e  uv  u7  v  ' 

The  effective  plastic  strain  e  is  calculated  in  accordance  with  the  first  part  of 
equation  (Alb).  The  true  value  of  ct  ,  namely  a  ,  is  used  because  plastic  flow  in  the 
true  stress  -  conventional  strain  state  is  to  be  calculated.  Next,  the  effective  total 
strain  e  is  obtained  by  using  the  second  part  of  equation  (Alb)  with  the  constants  v 
and  E  evaluated  as  0.  333  and  10^  pounds  per  square  inch,  respectively. 

The  elastic  component  of  the  circumferential  strain  is  assumed  to  obey  Hooke’s 
law,  which,  for  a  membrane  sphere,  is 


e 


9e 


1  -  v 
E 


Q  9 


=  0.  667 


(C2) 


The  plastic  component  of  the  circumferential  strain  e^p  is  calculated  by  equations  (Ala) 
and  (A2b).  Since  =  e ^  in  a  membrane  sphere 


etV'o-V 


which  leads  to 


e  -€ep 
e0P  9 


(C3) 


The  total  circumferential  strain  e  &  is  equal  to  the  sum  of  the  elastic  and  plastic  com* 
ponents. 

The  transverse  strain  is  calculated  from  equation  (A3).  For  a  membrane 
sphere,  equation  (A3)  becomes 


or,  for  v  =  1/3, 


=  -(e 


9 e  +  e  ep^ 


(C4) 
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The  instantaneous  values  of  the  radius  and  wall  thickness  in  a  membrane  region  of 
the  sphere  are  calculated  as  follows: 

a  =  aQ(l  +eQ)  =  6.  03(1  +  e  Q)  (C5) 

and 

h  =  hQ(l +e^)  =  0.060(1  +e)  (C6) 

Finally,  the  internal  pressures  that  produce  the  calculated  stress  and  strain  states  are 
calculated  by  the  equilibrium  equation 

"ah 

P  =  2-f-  (Cl) 
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